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Abstract 

In this paper the old problem of determining the discrete spectrum 
of a multi-particle Hamiltonian is reconsidered. The aim is to bring a 
fermionic Hamiltonian for arbitrary numbers N of particles by analytical 
means into a shape such that modern numerical methods can success- 
fully be applied. For this purpose the Cook-Schroeck Formalism is taken 
as starting point. This includes the use of the occupation number rep- 
resentation. It is shown that the N-particle Hamiltonian is determined 
in a canonical way by a fictional 2-particle Hamiltonian. A special ap- 
proximation of this 2-particle operator delivers an approximation of the 
N-particle Hamiltonian, which is the orthogonal sum of finite dimensional 
operators. A complete classification of the matrices of these operators is 
given. Finally the method presented here is formulated as a work pro- 
gram for practical applications. The connection with other methods for 
solving the same problem is discussed. 



1 Introduction 



One of the central problems of many-particle quantum mechanics, if not its 
main problem, is calculating the spectral representation of a many-particle 
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Hamiltonian, which typically has the form 



Here Kj contains the kinetic energy of particle j and the external fields 
acting upon j, and Wjk is the interaction of the particles j and k. As is 
well-known, this problem has a solution if Wjk = 0. On the other hand, if 
Wjk does not vanish, the problem is "almost" unsolvable in a strict sense. 
But the situation is not hopeless. For, what is really needed for practical 
purposes, is a "good" approximate solution. 

In this last field a tremendous work has been done, both analytically 
and numerically. Its mainstreams are well-known under the labels Thomas- 
Fermi method ([T],[2]), Hartree-Fock method ([3], [4]), density functional 
theory ([5|,[6|), configuration interaction method, Haken's method and oth- 
ers. With respect to these methods and their applications and refinements I 
refer e.g. to the following books [7], [8], [9j. There in addition an abundance 
of papers and monographs is cited, where the methods are also described in 
detail. 

A common feature of these procedures is that they contain one step in 
which a one-particle approximation of the N-particle problem is carried 
through. With the methods of Thomas-Fermi and of Hartree-Fock it is all, 
what is done. With the other methods the described first step is followed by 
other ones thereby improving the accuracy of approximation. Especially, by 
combining analytical and numerical mathematics great progress is achieved. 
Today problems can be solved which were regarded as unsolvable few decades 
ago. 

Nevertheless, the question is obvious, whether there are other approaches 
to a solution of the A'^-particle problem in quantum mechanics than those 
mentioned above. It is the aim of this paper to present such a new procedure. 
For this purpose I need some mathematical tools which, though they are 
widely known, I have briefly described in Appendix A.l. In particular the 
reader will find all the notation which is used throughout the text. (More 
details can be found in [lOj . The basic idea of the procedure as well 

as the main results are sketched in Section 2.3. 

2 The Structure of A^-Particle Hamiltonians 

2.1: In what follows only systems of particles of the same kind are consid- 
ered. When one starts studying a concrete sytem, its Hamiltonian is usually 
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defined using the position-spin representation, i.e. the Hamiltonian is an 
operator in the Hilbert space 

(g)^(L2(M3) (g) where (g) and (g) denote 
tensor products, and where is the complex vector space of spin func- 
tions (cf. Section A. 2.1). For exphcit calculations this representation is very 
useful. But the aim of this paper is primarily a structural analysis of the 
Hamiltonians of a certain class of systems, and in this case a more abstract 
formalism is adequate. It turns out that the Cook-Schroeck formalism (cf. 
Appendix A.l) is very useful for this purpose. 

Then our starting point is an arbitrary initial Hamiltonian of the shape 
(1.1), which is denoted Hn and defined in a Hilbert space := H^, 
where is the Hilbert space of the corresponding one-particle system. 
Now let K be the operator defined in "^^^ which contains the kinetic energy of 
one particle of a certain kind together with the action of the external fields. 
Moreover, let W be that operator in which represents the interaction 
of two particles of the kind considered. Then, using Formula (|A.1.25|) . Hn 
defined in is given by 

(2.1) HN = nN{K) + nN{w), 

where 

n^iK) := {{N - !)!)-! Y.PeS^ U{P){K l)U\P), 

(2.2) 

^n{W) := (2(iV - 2)!)-i Epg5^ U{P){W 1 ... l)U%P), 

and U{P) is the unitary permutation operator defined by the particle per- 
mutation P. Thus, using Formula (IA.1.27p . the operator i/^r specified for 
Bosons or Fermions reads 

(2.3) H^ = n%{K) + ^%(W). 

Here the definition := S^AS^ for an arbitrary operator A in is 
applied, where is the symmetrizer (+) resp. the antisymmetrizer (-). 
Then A^ is defined in the Hilbert space Ti^ = S^[W]. 
It is well-known that the structure of Hj^ given by (j2.3p is not helpful for 
studying its spectral problem, because the operators i}^{K) and i}^{W) do 
not commute. This suggests the question if it is possible to find an operator 
T acting in n'^^ , 1<M<N such that 

(2.4) = n^{T). 

Because the two-particle operator W cannot be represented by a one-particle 
operator it holds that M > 2. If the Hamiltonians as well as the operators 
K and W are selfadjoint, it turns out that M = 2 is possible as shown by 
the following 
Proposition 2.1: Let 



(2.5) 



H2h) = l{K ®l) + -i{l®K) + W 
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so that -^2(7) is defined in H'^, and let 70 := {N — 1)^^. Then 
(2.6) H^ = n%{H2ilo)) and H^ + ^^^(^2(7)), 7 + 70- 

Proof: Using (IA.1.28j) yields 

^%{K) =NS%{K®\® ...®\)S% 



(2.7) 



iV- 1 V 2 



^N~i (^) ^^^^^ ® i^) ^5 . . . ® 1)5± 



1 



Since by supposition H2{'y),K and are selfadjoint, ()A.1.30P can be applied 
so that with the help of (j2.7p the following relation holds: 



f]± (^2(7)) D 7f^^(i^ ^5 1) + 7^^^(1 + n%{w) 

(2.8) 

= 7(iV-l)f^^(i^) + J^^(H^). 



The term in the last line of (j2.8p is selfadjoint because it is the Hamiltonian 
of a (possibly fictional) A^-particle system. Since also 0^(^2(7)) is selfad- 
joint (cf. Proposition A. 1.11), relation ()2.8p is an equation, from which the 
proposition follows immediately. 

2.2: This result is somewhat surprising. The initial Hamiltonian H^ is not 
determined by H2{1), i.e. by a Hamiltonian of a system of two particles 
of the same kind, which is described by H^. Rather H^ is determined by 
-^2(70)) 7(7^ = — 1, which is a two-particle Hamiltonian for particles of 
mass (A^ — l)mo and external fields weakened by a factor (A^ — 1)~^, but 
with the same interaction W as the particles described by H^, which are 
supposed to have mass niQ. 

The system described by -^2(70) is fictional. I call it dummy system and the 
operator -^2(70) dummy Hamiltonian. 

In Appendix A2 two simple examples are given describing dummy helium 
and a solid with two dummy electrons. 

In what follows, the operator i72(7),7 7^ 7o is not needed anymore. There- 
fore it is convenient to use the notation i?2(7o) = -^20- 
Corollary 2.2: Because of (IA.1.3ip it fohows that 

(2.9) H^ = n%{H2o) = n%{Hi,) . 

2.3: Formula 1^2. 9\l suggests the basic idea of this paper: find an approxima- 
tion of H^ via an approximation of H^ such that the spectral problem of 
H^ can he solved approximately. 
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The details of this program are carried out for fermions in four steps, which 
correspond to the Chapters 3 to 6. 

Chapter 3 contains a formal analysis of a restriction of the initial Hamil- 
tonian . The aim is expressing the matrix elements of in terms of the 
matrix elements of the restricted dummy operator H2Q, which is bounded. 
The results are summarized in Proposition 3.10. 

In order to use them for the present purposes, a properly chosen orthonor- 
mal system Oi in the one-particle Hilbert space is needed. In Chapter 
4 arguments are given that such a system is obtained via the Hartree-Fock 
procedure applied to 5"^. Thus the restrictions used in Chapter 3 can be 
justified. Moreover, a heuristic argumentation suggests that can be 
"truncated" such that an operator results, which is, depending on a 
parameter a S N, an approximation of H2q- 

In Chapter 5 it is shown that the operators converge strongly to H2Q, if 
a — >• 00. This has the consequence that the operators = also 
converge strongly to H^. Therefore it is possible to apply the results of the 
theory of spectral approximation (cf. e.g. |12j. [13j ). 

Finally, in Chapter 6 an analysis of the operators is given. It is shown 
that they are block-diagonal, i.e. their matrices with respect to the cho- 
sen orthogonal basis are orthogonal sums of finite dimensional matrices, the 
structure of which are analyzed in detail. At this point numerical methods 
can come into play. 

In Chapter 7 the results of the previous chapters are summarized in the form 
of a work program, which can be regarded as the main result of this paper. 

3 The Hamiltonian Hj^ and its Matrix 

3.1 Preliminary remarks 

3.1.1: Since in what follows only such systems are considered which consist 
of fermions of the same kind, the notation introduced in Appendix A.l can 
be used throughout. The Hamiltonians of these systems usually have the 
following property. 

1. ) They are unbounded, but bounded from below. 

2. ) Their spectrum below a certain value eo is discrete, otherwise continuous 
with possibly inserted discrete values. 

Since in this paper we are only interested in the discrete spectrum outside 
the continuum, i.e. in the bound states of the system, the question is ob- 
vious, whether it is possible to restrict the spectral problem to the discrete 
eigenvalues. In other words: is it possible to realize the following 
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Assumption 3.1: There is a subspace Ti^ C Tl^ such that the restriction 
H2Q of H20 to the subspace T-L^ C H'^ is bounded (so that it can be defined 



'20 -'-'20 

on 'H'^), and has the same discrete eigenvalues as outside its continuous 
spectrum. 

In this chapter it is assumed that such a subspace Ti^ C exists. Then 
H^ = r2^(i/^) is bounded and defined on T-L^. This operator is the subject 
studied in the following sections. In chapter 4 arguments are given that the 
assumption can be realized. 

3.1.2: The starting point for the further considerations is the following 
Notation 3.2: 1.) Let ^1 = {(/j^ : ^ G N} be an arbitrary ONB in H^, and 

let 

(3.1) (t>.,....M := (^.i ^ • • • ® G 

with 2 < M < and G N,j = 1, ....,M. Then an ONB C Ti^^ is 
defined by the vectors 

with S^,j being the antisymmetrizer (cf. (|A.1.6p ). 

2. )For each sequence ki • • • hm of indices there is an infinite sequence k := 
{ki, k2, ks, ■ ■ ■ , ) of so called ocupation numbers k^ defined by 

M 

(3.3) k^ = '^6^,f,^. 

Hence /c^ = 1 or 0. Moreover there is a one-to-one correspondence 

(3.4) k i — > Ki ■ ■ ■ Km- 
Thus we can write 

(3.5) *± =: ^Uk)- 

3. )The term "sequence of occupation numbers" is abbreviated by bzf and 
the set of all bzf, which have exactly M numbers 1 is denoted BZFm- The 
set BZF comprises all hzf . 

Now let ( •, • )2 be the inner product in and let 

(3.6) E{k,m) := (^^(fc), 7^20^2 (m))2 

be the matrix elements of the dummy Hamiltonian. Then the matrix rep- 
resentation of H2Q reads: 

(3.7) i/2o = ^^i?(^,m)*2(fc)(M/2-(m), • )2 . 

k rri- 

Using the abbreviation 

(3.8) vI,2(fc)(M/-(m), •)2=:r(^,m) 
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together with Formula (|A.1.28p yields: 



(3.9) 



N 



( 2 ) Efc ECk, m)S],{TCk, m)0l0---^l)S^. 



Remark 3.3: Here and in what follows the sums and Ylm under- 
stood to run over all bzf which occur in the elements of Each of these 
sums can be arbitrarily ordered because each ordering of the k or the rh 
yields an ONB. Since is assumed to be bounded the sums J2k 
can be interchanged. 

Notation 3.4: 1.) As usual the abbreviation 

(3.10) (^'li^^vl") ■■= {^N{n'),H]:,^]^ih)) 

is used. 

2.) Let n G BZFj. In the present case j = 2 or j = N, and in the 
next section also j = M is used with 2 < M < N. But irrespective of 
these special choices, for each two bzf an addition and a subtraction can 
be defined by adding, respectively by subtracting their components. Since 
these operations on two bzf not necessarily result in a bzf the following 
notation is used (cf. IA.1.3P : "h ±1% is a bzf" or "n it m G BZF" . These 
expressions indicate that the sequence hzizrh does not contain the numbers 
2 or -1. 



3.2 The basic lemma 



3.2.1: It will be shown that the following proposition holds. 

Lemma 3.5: There is a function C such that for each triple (n, k, rh) of 

bzf with Yl,a ^a = ^,Yli3^l3 — J2(3 — ^ relations 

(3 11) C{h,k,7fi) =±1, if h — rh ^ BZF and h + k — rh £ BZFi\f 
= 0, if n — rh ^ BZF or h + k — rh ^ BZF]\j 

hold, and that moreover 

(3.12) {n\H^\n) = C{n, k, rh)E{k, rh)5{n , h + k — rh), 

k m 

where 6 is the Kronecker symbol and where E{k, rh) is defined by ()3.6|) (cf. 
also Remark 3.3). 

From Lemma 3.5 one can draw the following 

Conclusion 3.6: If h and h' are given, {n'\H^\h) can be unequal zero only 
if n — m G BZFm-2 and n' — k G BZFn-2. These relations can be satisfied 
only for {^) bzf m and k. Hence the sums in (j3.12p have finitely many 
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summands. 

3.2.2: Though the lemma is used in this paper solely in the above version, 
for later purposes a generalization of it will be proved in the next sections. 
(The expenditure is the same in both cases.) In order to do so, some notation 
is introduced. 

Let 2 < M < N and let Am be a bounded operator defined on T-LJ^ such 
that 

(3.13) An = n]^iAM) 

is defined on Ti^. The matrix elements of Am with respect of are again 
denoted E{k, m) so that here 

(3.14) J];A:, = ^m„ = M. 

a a 

Moreover let 

(3.15) {h'\AN\n) := {^]^ih'), AN^]^{h)) . 

Finally, the function C is defined as in Lemma 3.2 but with condition (|3.14p . 
Proposition 3.7: The relation 

(3.16) {n\A]\[\h) = C{h, k, m)E{k, m)6{h' ,h + k — rh) 

k rh 

holds. (Cf. also Remark 3.3 and Conclusion 3.6.) 

3.3 Proof of Formula (3.11) 

3.3.1: The starting point is Formula (1A.1.28P and the analogue to Formula 
^Mi- Thus 

(3.17) {n'\AN\n) = f ^ J ^ J] E{k, rh)Z{n , n + k-rh), 

k m 

where 

(3.18) Z(n', n, k, rh) = (*^(n'), (T(fc, m) ® 1 • • • ® l)^'^(n)). 

Here the operator T(k,m) is defined by strict analogy with ()3.8p . Thus 
(3.19) 

T{k,m) = ^-,(fc)(M'-,(m), ■)m 

= (M!)-lEQe5„^"(Q)^A/(^)(</'Mp-ia)-A'«-i(M,' • )m , 
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where fJ-i, - ■ ■ , /iAf with /^i < . . . < /iAf is the correspondence defined by 
(|A.1.13p . Now using the correspondence ui, . . .uj^ on with h'l < . . . < 
one finds that 

(3.20) {Tik, m) 1 ... ® l)^'^(n) = ^'^^(A:) x{n, m), 

where 
(3.21) 

x{h,m) = (mM\)-^2 a-{P)a-{Q) 

P&Sn 
QeSM 

'^^i=l^'^/^<3-iO)''^^p-i(j))l)'^'^p-i(M+i)-^p-i(jv)- 
Thus finally we obtain the relation 

(3.22) Zih', h, k, m) = (^]^(n'), ^Ijik) ^ x{n, m)). 
3.3.2: In this subsection the following proposition is proved: 

(3.23) x{n,m)^0, 
if and only if n — m G BZF. 

Firstly it is assumed that h — rh ^ BZF. Then there is a number a such 
that rria = 1 and = 0. Consequently, for each Q G Sm there is an r for 
which Q~^(r) = a holds. But for each P G Sn the relation P~^(r) 7^ a is 
true. Thus for each pair P, Q 

so that also x(^) "^) = 0. 

Secondly let us assume that fi — rh£ BZF. Then for each a with rria = 1 
also Ua = 1 holds. Consequently one has to look for all pairs Q, P such that 

For all other pairs Q, P the product in (|3.25p is zero because {(j)^, 4>u)i = S^iu- 
Thus (j3.25p is equivalent to 

(3.26) ^Q-i(i) = ^P-Hi)^ J = Ir • • ,M. 

In order to satisfy (j3.26p . for a given Q G Sm the permutation P G Sn must 
be such that the /xi, • • • , ^m, which by presumption occur in z/i, • • • , vm^ 
occupy the places 1, • • • , M being ordered by Q. All these pairs Q, P can 
be explicitly indicated by the following procedure. 

Let S G Sn be that permutation for which 

(3.27) (i^s'-i(i)' ■ ■ ■ --^s-^iN)) = ifJ'ii--- ■ifJ'M-.Qi-:- ■ ■ ,Qn-m), 
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where Qi, - • • , Qn-m are all those i^i, • • • , i^at which are unequal fii, • • • , fiM, 
and in addition let < • • • < qn-m- Then with the help of (|A.1.7p one 
finds 

^^(n) =a-{S)^/N'\S-^U{S)(p,,...,^ 

(3.28) 

= cr-(5)\/iV!S'j^</'/.i-/.Mei-eiv-M 

so that (j3.2ip now reads 
(3.29) 

X(n,m) = {N\M^r^^a-{S)Y.p,Q<y~{P)<y-{Q) 

It follows from (I3.29P that only those pairs Q, P give nonzero summands for 
which P has the form: 

(3-30) P=(q, Z + \n ) ( I,--. |m 



where R G Sn-m is an arbitrary permutation. Hence, for a given Q there 
are {N — M) \ permutations P of the form (I3.30p such that ()3.26p is satisfied. 

Since R acts on (i^5-i(m+i), • • • i^5-i(iv)) = iei,--- ,Qn-m) one finally 
obtains 
(3.31) 



X(n» = {N\M\)-2a-{S)Y.^Qa-{R)a-{Qfct>,^_,^^^...,^_, 



(N-M) 



for all pairs n, rh with h — m & BZF. Hence the proof of relation (j3.23p is 
complete, and in addition the explicit form of x(n, rh) is obtained. 
3.3.3: Inserting (fS^TD into (15:20]) and (l3T8]) yields 
(3.32) 

/N\^^ _ _ _ _ - 

Z{h',h,k,m) = l\ a {S){^^{h'),S^{^j^j{k)(g)^N~Min-m))). 

With the help of ()A.1.7p and by the correspondence k o {ki, • • • ,km) it 
follows that 

(3.33) = (M!(iV - M)\)-^S]^{<P,,...,,,, ^ K-bn-m) 

= Q-K-iT)^-^{h + k-m), 
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where T £ Sjy is the permutation which hnes up the sequence • • • , km, 
£"!)■■■ ) QN-m) in its natural order. Thus one obtains 

Z{h', n, k, m) = 
(3.34) = (f,)"V-(5 • T)(^^(n'), ^^(n + k-m)) 

= {My^<y'{S-T)8{h',n + k-m). 

It follows that Z{n' ,h, k,m) ^ exactly if n — m G BZF,h + k — rh G EZFj^ 
and n' = h + k — m. 

3.3.4: Since the permutations S and T are uniquely defined by the sequences 
of indices (z^i, • • • , I'm)-, (/^i, • • • , /^m) and (ki, . . . , km) or equivalently by 
n, m and it is obvious to define the function C by 
(3.35) 

C(n, /c, m) = a~{T ■ S) = ±1, if fi — rh € BZF and fi + k — m £ BZF]\f 
= 0, otherwise. 

Now inserting (I3.34p together with (I3.35|) into (I3.17P Formula (I3.16|) is seen 
to hold, thus also Lemma (3.5). With respect to the sums and Yl,m ^ 
refer to Remark 3.3 and Conclusion 3.6. 

3.4 An algorithm for C{h,k,rh) 

3.4.1: The question to be answered in this section reads: is there a finite 
procedure for calculating C(n, A;,m) if h,k,m are given bzf. As in Section 
3.3 the more general case 2 < M < A is considered. 

Since by definition C(n, k, rh) = if the condition 

(3.36) n — m £ BZF and h + k — rh G BZF]^ 

does not hold, only the case needs to be considered that ()3.36p is true. Then 

(3.37) C{n,k,rh) ■.= a-{T-S) = (-l)'^(^)(-l)-^(^), 

where S is defined by (|3.27p and T by (j3.33p . Moreover, J{P) here means 
the number of inversions of a permutation P (cf. e.g. (jA.l.lOP ). 
3.4.2: To begin with, J{S) is to be calculated. Let h be given. Then 
exactly A numbers = 1, • • • , A exist such that rii^- = 1 and ui < ■ ■ ■ < 
UN- Hence n o {vi ■ ■ ■ ^v^)- Likewise, if m is given, exactly M numbers 
//j, j = 1, • • • M exist such that m^. = 1 and /ii < • • • < fiM- 

Because of n — m G BZF for each j G {!,••• M} there is an rj such that 

(3.38) = o-nd j < rj . 



12 A NEW APPROACH TO THE N-PARTICLE PROBLEM IN QM 



The permutation S is defined by (|3.27|) . i.e. 

(3.39) (l^5-i(i)), • • • , l^s-i(Ar)) = (/^l, • • • ,/iM, Ql,--- , QN~m) 

and Ql < ■ ■ ■ < qn-m- Then the right-hand side of (|3.39p can be generated 
from (z^i, • • • , vn) by the following procedure. 

First, /ii = Vj.^ is positioned at the ri — th place in (z/i, • • • , vn)- Therefore 
one needs ri — 1 inversions to bring /^i at the first place. Thereby the 
positions of ;U2, • • • , /^m in (z^i, • " " ■, i^n) are not changed. 

Second, ^2 = z^r2 is positioned at the r2*^ place in {vi - ■ ■ , vn) so that one 
needs r2 — 2 inversions to bring ^2 at the second place. Again the positions 
of /U3, • • • , /ijvf are not changed. 

Thus, in order to bring /ij = u^j to position j one needs rj — j inversions. 
Therefore the total number of inversions, which realize the permutation S 
in (]3.39p . is given by 

M 

(3.40) j(5) = ^^.__M(M + l). 

i=i 

3.4.3: Now, J{T) is to be determined. This task is the following. Let h — rh 
and k be given. Then h — rh ^ BZF corresponds to the sequence of indices 
(^^i)""" tQn-m) and k to the sequence (ki,--- ,km,)- Thus the sequence 
of indices (ki, • • • , ka/ , £»! • • • , Qn-m), brought to its natural order by the 
permutation T and denoted {i^'i, - • • , i^^), corresponds to n + fc — m £ BZF. 
Hence 

(3.41) (i^T(i))' ■ ■ ■ ' '^r(Af)) = (/^i) ■ ■ ■ ) i^M, Ql,-- - , Qn^m)- 

Therefore, because n + k — rh is given, also (i^(, • • • ^v'^) is determined so 
that for each j G {1, . . . , M} the position Sj of in (i/[, • • • , v'^) can be 
read off with the help of the relations 

(3.42) Kj = v'^^ and j < sj . 

Using the same arguments as in Section 3.4.2 one obtains for T~^ the result 

AI 

(3.43) J(r) = J(T-i) = ^s,- -M{M + 1) . 

i=i 

3.4.4: Finally, the algorithm for C(n, k ,m) can be formulated thus: 
1** step: Take bzf h, k and rh which fulfil the equations = 
Y^ka = J2a^a = and test Condition (j3.36|) . If it is satisfied go to 
the next step. If it is not, define C(n, k,7n) = 0, so that the task has been 
done. 
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2 step: Take n, m and determine the corresponding sequences of in- 
dices (ui, - ■ ■ , un) and (//i, • • • , hm)- Then from (|3.38|) read off the numbers 
rj,j = !,••• , M, and calculate J{S) with the help of ()3.40p . 
3^^^ step: Take h + k — m and k, and determine the corresponding se- 
quences - • • ,i^'j\f) and • • • , km)- Then from (I3.42|) read off the num- 
bers Sj, J = 1, • • • , M and calculate J{T) with the help of (|3.43|) . 
4*^* step: Calculate C{h,k,m) using (|3^ . 

The coefficients C{h,k,rn) do not depend on the specific physical system, 
for which they are used, rather they are completely combinatorial. In other 
words, they result solely from the algebraic structure imposed on the set 
BZF. Therefore they can be computationally calculated once for all. A 
trivial special result is the following. 
If /c = m, then T^^ = S so that 

(3.44) C{h,m,m) = 1. 

3.5 The final form of {h'\H^\h) 

3.5.1: For the sake of simplicity in this section only the special case M = 2 
is considered. This does not entail any loss, because the Hamiltonians we 
are interested in this paper are supposed to have two-particle interactions. 
The starting point for this section therefore is (I3.12p . Moreover it is assumed 
that the matrix elements E{k, rh) of the dummy Hamiltonian Hi^Q and the 
coefficients C{n,k,m) are given. 

Now the problem to be solved reads as follows. Let the pair h' , h be given. 
Then determine those pairs k, rh for which the summands in and 
do not vanish on general grounds. 

It is known from the previous considerations that for given n', n only those 
k and rh in (I3.12P are relevant which satisfy the condition 

(3.45) n — TTi G BZF, h + k — rh G EZF^, h' = h + k — rh, 
hence also 

(3.46) n' -k € BZF, fi - h = k - rh. 

Consequently, the sums in (|3.12p have only finitely many summands as al- 
ready remarked in Conclusion 3.3. 

3.5.2: In this section a disjoint dissection of all possible pairs k, rh for given 
h' , h will be defined. For this purpose it is useful to introduce some new 
Notation 3.8: Let k,rh € BZF2 so that Yl^a ~ = 2 is satisfied. 
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Then the sequence d := k — rh is cahed differences sequence. The set of all 
difference sequences is denoted V. 

Consequently there are only three types Vg, g = 0,1,2 oi possible d gen- 
erated hy d = k — m : 

T>Q contains only one element 6 := (0, • • • , 0, • • • ), and 6 is generated by all 
k, rh with k = rh. 

T>i contains d = (1, —1, 0, • • • ) and all permutations of it. They are gener- 
ated by all k,rh such that for exactly one a the relation k^ = rua = 1 holds. 
P2 contains d = (1, 1, —1, —1,0, • • • ) and all permutations of it. They are 
generated by all k, rh for which no a exists such that ka = rua = 1. 
It follows immediately that the sets 'Dq,'Di,V2, are disjoint and that 

(3.47) P = Po,UPiUp2- 

Many results of the next sections and chapters are based on the following 
Proposition 3.9: For each pair n', h G BZF^ there is a set 

(3.48) {di,...,dL] dVi 
such that 

(3.49) n' = h + di + ... di- 

The number L is uniquely determined by h' and n, but the difference se- 
quences dj, j = 1, . . . , L are not. 
Proof: From the pair h', h one forms the matrix 

(3.50) X:=( t Uf^^' 
For each column of X the following alternative holds: 

(::) ^ Q - " < 

Because n and h' both contain N numbers 1, there are equally many columns 
Cj*) and (q) . Let L by the number of each of the two kinds. Moreover let 
Qj and aj, j = 1, - ■ ■ ,L he numberings of the indices of these columns such 
that 

,3.52) (;;,.).(;) and K^ 



Then for each pair of indices Qj,(Tj,j = 1, • • • ,L define dj = {dji,dj2 
by djg. = l,dja = - 1 and dja = 0,a / Qj,(Tj. 
It follows from (13321) that 

(3.53) h + dj = {■■■ , n„, • • • ,n'^^,--- ,n'„^,--- ,np,---) , 
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if Qj < Gj , and analogously, if CTj < gj . 

Since j 7^ i implies gj / gi, aj 7^ ai the addition of di to n+dj can be carried 
through without altering n'^. and n'^. in (j3.53p . Thus, finally one ends up 
with (|3.49p so that the proposition is proved. 
Three immediate consequences are useful later on. 

1. ) For each pair z, j with i ^ j the relation di + dj G P2 holds. 

2. ) h' — h £ Vg, exactly if L = ^j, g = 0,1, 2. 

3. ) h' — fi ^ T> exactly if L > 3. 

3.5.3: Using the results of the previous sections the problem formulated 
in 3.5.1 now can be solved by giving a disjoint classification of the matrix 
elements defined by (|3.12p . According to (|3.47p four cases have to be taken 
into account. 

1*** case: h' — h ^ D. It follows immediately from (I3.12p that 

(3.54) (n |F^|n) = 0. 

2nd case: h' — h £ Vq, i.e. h' = h. Then {n\Hjj\n) is unequal zero only, 
if also k = m. Hence the double sum Ylm I'educes to a single sum Ylm.- 
This sum runs over all rh for which h — rh £ BZF holds. Since h contains 
N numbers 1 there are exactly (^) different sequences rh such that this 
condition is satisfied. Because of C{h,rh,rh) = 1, define 

(3.55) 0) = E{rh, rh) 

rh 

for all rh with h — rh £ BZF. Thus finally 

(3.56) {n\H^\n) = £{n,d) . 

3^^^ case: h' — h = di £ Vi. Hence h + di £ BZF. Then {h'\H'^\h) is 

unequal zero only, \i k = rh + di. Therefore the double sum Ylm ^-g^iii 
reduces to a single sum which runs over all rh so that h — rh £ BZF 
and rh + di £ BZF. These rh can be characterized as follows. 

Let di be given by di^ = l,di^ = — 1 and dijs = 0, /? 7^ k, /u. Hence 72^ = 
and = 1. Then rh + di £ BZF if and only if = and = 1. In 
order to satisfy the condition h — rh £ BZF it is necessary and sufficient 
that = 1 and that there is an a 7^ k, ^u, for which ria = 1 and = 1 
holds. Since = and = 1 there are — 1 numbers a 7^ k, ^ for which 
= 1 so that runs over all rh for which rria = = 1. Now define 



(3.57) <?(ri, di) = C{'h,rh + di,rh)E{rh + di,rh). 

rh 
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Then 

(3.58) {h'\H]^\h) = £{n,di), 
if n' = n + di , di £ Vi 

4*^ case: h' — h = d2 G 'P2- Let ^2 be defined by ^2^ = d2A = li d2fi = d2u = 
— 1 . Then k and rh with k — 7fi = d2 are uniquely determined hy k^^ = k\ = 1 
and Tfi^ — Tfiy — 1. Hence the ^^duces to a single term. Now define 

(3.59) £{n, d2) = C{h,rh + d2,rh)E{7h + d2,m) , 

where m is determined by = = 1. Then, if h' — n = d2 £ T>2 

(3.60) {h'\H-\n) = £{h,d2). 
3.5.4: Summing up one obtains 

Proposition 3.10: 1.) The matrix elements of the fermionic Hamiltonian 
H]^ defined in Assumption 3.1 (Section 3.1.1.) are given by 

. , {h'\H^\n) =£{h,h' — h), if h' — h £ V 

=0, if h'-n^V 

and 

(3.62) £{h\h — n) = £{h,n — h). 

2.) Let be given fi £ BZF and d £ D so that h + d £ BZF. Then it follows 
from Formula ()3.12p that 

(3.63) iS(n, d) = C(n, rfi + d, m)E{m + d,m) , 

m 

where the sum runs over all m for which n — rh £ BZF and m + d =: k £ 
BZF2. Hence, the matrix elements {h'\HJj\n) are determined solely by the 
matrix elements E{k, rh) of -f^20' which k — rh = d = h' — h. The sum in 
dm]) is finite. 

For bosons a result holds, which is formally equal to (|3.12p . (|3.6ip and 
p.63p . but the terms C,E,£ are defined differently. 



4 Heuristic Considerations 



4.1: The result of Chapter 3, summarized in Formula (j3.6ip . is completely 
formal up to now. This is due to two unsolved questions connected with it. 
They read as follows. 

First, can Assumption 3.1 he verified? More concretely, is it possible to find 
a one-particle Hilbert space T-L^ such that the dummy Hamiltonian H2Q is 
defined on H.^ and is bounded, and such that in addition the spectrum of 
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contains the discrete eigenvalues of the dummy Hamiltonian H2Q, which 
is primarily defined ? (Cf. e.g. (|A.2.ip . (|A.2.4p . (|A.2.6p .) 
Second, is Formula 13. 61\) of any advantage for the spectral problem of ? 
4.2: To begin with, let us look for an answer to the first question. As is 
already described in Subsection 2.f, any investigation of the Hamiltonian 
Hn of an A'^-particle system starts with a more or less informal specification 
of the external fields acting upon the particles and of their interactions. 
Customarily this is done using the position-spin representation. Then Hi\i 
is of the form (jl.ip or, what is the same, (j2.ip . It is densely defined in a 
Hilbert space . Likewise the corresponding dummy Hamiltonian can 
be immediately written down as is shown for two examples in appendix A2. 
It is defined in a dense linear submanifold of H^. 

In order to verify Assumption 3.1 the space has to be properly re- 
stricted to a subspace T-l^. Such a restriction in turn can be carried through 
by finding a proper orthonormal system Oi in so that = span Oi. 
Having in mind the physical meaning oi TH} suggests taking the Hartree- 
Fock procedure for .ff^ to determine Oi. This is because this procedure is 
based on the Ritz variational principle which guarantees optimal approxi- 
mation. Disregarding the fact that the Hartree-Fock procedure generally is 
an infinite task, let us assume that it is completely carried through for the 
dummy Hamiltonian in T-L"^. Thus one has obtained an orthonormal 
system O2 C 'H'^ of vectors 

(4.1) = - (Aa ^ K < A 

together with corresponding energy levels E^^x = E(m, m) for k, A m, 
which approximate the discrete eigenvalues of H2Q outside its continuous 
spectrum. Since the Hamiltonians considered in this paper are supposed to 
have a bounded discrete spectrum, the set of the energy levels E^\ is also 
bounded. 

Now define Oi to be the set of all (j)^ for which there is a ^'^^ in 
Hence V} C V} is defined and d is an ONE of V} . 

Then the restriction //go of ^20 *o space 'H^_ is bounded because its 
spectrum is approximated by the set {-E^A : k < A} and its eigenvectors by 
the set 02- As usual can be defined on the whole space using its 
matrix representation with respect to the ONB 02- 

In most cases of practical application the complete Hartree-Fock procedure 
cannot be achieved, because it is infinite. Therefore one has to content one- 
self with a finite section of this procedure. But also such a finite procedure 
can be complicated. 
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Thus other methods were invented which are equivalent to the Hartree-Fock 
procedure or approximate it. 

4.3: Let us now come to the second question. From (j3.6ip one draws 
immediately the simple consequence. 

Proposition 4.1: The matrix of the Hamiltonian H^ = Q^{H2q) is diag- 
onal if the matrix of the dummy Hamiltonian H2Q is diagonal. 

Unfortunately this result cannot be used to obtain the exact discrete 
eigenvalues of a realistic N-particle system, because the exact eigenvectors 
of a dummy Hamiltonian iif^ containing interaction are not elements of any 
ONB 82- But, if one contents with a Hartree-Fock approximation of the 
eigenvalues of H2Q, Proposition 4.1 delivers a Hartree-Fock-like approxima- 
tion for the eigenvalues of H^. 

Thus, if one wants to obtain better approximations, one has to solve the 
following problem. Since the results of Chapter 3 are valid for arbitrary or- 
thonormal bases Bi C Ti^, one firstly has to choose such an ONB and one 
has to calculate the matrix elements E{k,rh) of for the ONB 82- The 
second part of the problem then is the question, wether the choise of Bi is 
helpful for a reasonable approximation of and also of T^]^- 
A heuristic idea to cope with this problem is the following. Choose the 
ONBs Bi, resp. B2 such that the matrix E of with respect to B2 , 
i.e. the matrix defined by the elements E(k,m) is "as diagonal as possible". 
This means the elements of B2 should approximate the eigenvectors of 
optimally. Hence we end up again with the Hartree-Fock method, an equiv- 
alent of it or an approximation. Therefore Bi = Oi and B2 = O2 with Oi 
and O2 being defined as in Section 4.2. 

In order to get further insight into the general structure of the matrix 
E I will give some purely heuristic arguments, which are based on physical 
intuition. For this purpose let the Hartree-Fock energy levels E^\ be num- 
bered such that Efj_^ < E^\ if // < 1/ < A. Then, if A is large enough, i.e. 
if it exceeds a certain value a one expects that the particle having state A 
is "almost" free, so that the interaction between the two particles having 
the states /j, and A is "almost" zero. Thus the two particles with states /x 
and A are "almost" free, if one particle of this pair is "almost" free. This 
implies, that the vector ^ "almost" an eigenvector of -^20- Now 

let /X, A correspond to a sequence of occupation numbers rh. Then E{m, rh) 
is "almost" an eigenvalue of iif^ for the eigenvector = ^'^(m), so 
that E{k,m) is "almost" equal to E(m,rh)6{k,m.). Consequently E{k,ifi) 
is "small", i.e. it is "almost" zero, if /c 7^ m. 



JOACHIM SCHROTER 



19 



Now let us suppose that the term "small" has been concretized. Then 
the above considerations can be summarized in the following assumption. 

There is a natural number a such that E(k, m) is small for each pair k, rh 
with k ^ rh, for which a k^j = l,a > a exists or an rUg = \, g > a. 

This suggests truncating the matrix E by substituting zeros for its small 
elements. Then intuitively one conjectures that the truncated matrix leads 
to an approximate solution of the spectral problem of HJ^ which is the goal 
of this paper. The precise meaning of the conjecture will be given in Section 
5.1. 

4.4: The first step is defining the truncated dummy Hamiltonian and 
drawing some consequences. The operator is determined by its matrix 
E, which is given by the elements: 

E{k, m) =0, if A; 7^ m and if a, g exist such that 

(4.2) fco- = 1, cr > a or rUg = 1, g > a, 
E{k,rh) = E{k,rh), otherwise. 

According to this definition is a finite nondiagonal matrix of order (2) 
with an infinite diagonal tail of elements E{rh, rh) where rh contains at least 
one rUg = 1, g > a. 

Consequently, is bounded so that H]^ = Q]^{H2o) is also bounded 
and has the form (j3.6ip . but with S{h,d) defined by E(k,m) the same way 
as £{h,d) is determined by E(k,m). Then we obtain the following 
Proposition 4.2: Let us consider {h'\H^\h) . 

1. ) li h' — h ^ D, one concludes from (|3.12p that h' — h ^ k — rh which 
implies (cf. ^M)) 

(4.3) {h'\H]^\h) = {h'\H]^\h) = 0. 

2. ) If h' — h = d £ D and if there is A > a such that dx 7^ 0, it follows from 
d = k — rh that fc;^ = 1 or mx = 1. Hence E{k, m) = 0, and consequently 

(4.4) {n'\H]^\n) = 0. 

3. ) Now let h' — h = d G V and suppose that = 0, if A > a. Then 
if d £ the bzf k and m are uniquely determined and one obtains (cf. 
Formula (I339D ) 

(4.5) {h'\H-\h) = i{h, d) = £{h, d) = {h'\H-\h). 

If d G Pi, things are more complicated. One has to apply the general method 
described in Section 3.5.3 for several special cases. The matrix element 
{h'\Hjj\h) can be zero or unequal zero, and it need not be equal to {h'\Hjj\h). 
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4.5: By these considerations the problem of determining the eigenvalues of 

H^ is transformed into the following two ones. 

First, do the eigenvalues of H^ approximate those of H^ ? 

Second, can the eigenvalues of H^ he calculated, at least partially? 

The first problem is treated in Chapter 5, the second one in Chapter 6. 



5 Spectral Approximation 

5.1: The dummy operators and -ff^ are defined on T-i'^ and are bounded 
so that D := — is also bounded and defined on H^. Moreover 
and D depend on the fixed number a. In what follows these operators are 
understood to be functions of a parameter a G N and a > 2 so that they 
are written D{a) and -ff^(a), and for the matrix elements of E we write 
Ea{k,m). Then by definition 

(5.1) {k\D{a)\m) = E{k, m) - Ea{k, m), a>2. 

Therefore 

{k\D{a)\m) =0, if k = m or k -H- ki, K2 < a and m -H- ^1,^2 < 

{k\D{a)\rri) = {k\D{2)\rh), ki > a or K2 > a 

or /ii > a or ^2 > 

These properties have the following consequence. 

Proposition 5.1: The sequence {D{a) : a € N, a > 2) converges strongly 
to 0. 

Proof. Let the projections F^ and F'^ be defined by 

00 

(5.3) E •) 

and F^ = 1 — Fa . Then Fa converges strongly to 0, if a — )■ 00, and F^ 
to 1. By a simple calculation using Formula (j5.2p one verifies that for each 

f^nh 



(5.4) 



FaD{a)f = FaD{2)f, 
F'aD{a)f = F'aD{2)Faf. 



Therefore 

(5.5) ma)f\\ < \\FaDi2)f\\ + \\F^D{2)Faf\\. 
Because of 

(5.6) ||F^Z?(2)F„/||<||Z)(2)||||F,/|| 
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and because converges to 0, the proposition is seen to hold. 

Now, in order to transfer the last result to H^{a) := Qj^{H2o{a)) one 
needs the following theorem. 

Proposition 5.2: Let be given a sequence {AM{r) : r G N) of bounded 
operators defined on T-L^. If A]\j{r) converges strongly to Om for r — )• oo, 
then also 

(5.7) s - lim nj^iAMir)) = 0^ 

r— oo 

i.e. 0^(AAf(r)) converges strongly. 
Proof: 1.) Let g en^, then 

9= ^ K 



If 



it follows that 



and 



Kl...ftJVf 



|2 

tM+l---fSjV I 

KJV/+1---KJV 



If S is a bounded operator on 'H*'^, then 



(5.8) ||(i?^l»...$3l)9f = Yl II^W+i. 



||2 



2.) Now, let us consider the operators AMir). By supposition AA./(r) con- 
verges strongly to Om- Hence by the principle of uniform boundedness (cf . 
e.g. |13j . p. 150) there is a number K such that for all r G N: 

\\AM{r)\\<K. 



Thus, one obtains for all r G N: 

(5-9) Y Pa/(^)Xk f < J2 

KJV/+1---KJV Ki\f+l---K]V 

Hence, by the criterion of Weierstrafi the left series converges uniformly with 
respect to the variable r. Therefore the limit r — )• oo can be interchanged 
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with the sum so that 

hm^^oo \\{AM{r) 1 • • • l)^^ 

(5.10) 

= E^M+i-^N lim^^oo \\AMir)xKM+i...Kj\'^ = • 
Thus, because for all g G T-L^ the relation 

\\n],iAM{r))g\\ < \\{AMir) 1 . . . )<7|| 



M 



holds, Formula (j5.7p is proved. 



Then, with the help of the Propositions 5.1 and 5.2 one obtains the fol- 
lowing consequences. 

Proposition 5.3: 1.) The sequence {H^{a) : q G N, q > 2) converges 
strongly to H^. In other words, the operators H~^{a) approximate H~^ (cf. 
|12) . p. 228, Formula 15. ip . Hence all results concerning the approximation 
of one operator by a strongly converging sequence of operators are valid for 
H^ and H^{a). (Cf. |12j . Chapter 5.). Here only two of these properties 
are sketched. 

2. ) If e is an isolated eigenvalue of H]^, then there is a sequence (e^ : a G N) 
such that ea is an eigenvalue of H^{a) and such that Ca e. (Cf. |12] . p. 
239, Theorem 5.12). 

3. ) Let P be the spectral projection belonging to an eigenvalue e of H^. 
Then an uq exists such that for each a > ao there is a spectral projection P^ 
belonging to H^{a) and such that Pa converges strongly to P for a — )• oo. 
(Cf. [12J, p. 240, Theorem 5.13.) 

5.2: Besides the approximation of H^ by H^{a) sketched above there are 
results, which are based on the norm of D{a), i.e. on 6{a) := ||D(a)||. Using 
the Formulae ()A.1.29p and ()A.1.30p one obtains the relation 

(5.11) \\H--H-{a)\\ =11 n-{D{a)) \\< (^)<5(a). 

Thus it follows that H^{a) converges in norm to H'j^ if 

(5.12) lim 6{a) = 0. 



Theorems concerning spectral approximation based on convergence in norm 
can be found in [12], p. 291, Theorem 4.10, p. 362, Theorem 5.10 and in 
[13] . p. 249, Proposition 5.28. 
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6 The operator H]^ and its matrix 
6.1 Preliminary remarks 

6.1.1: The aim of this chapter is proving 

Proposition 6.1: 1.) The operator defined on T-L^ is an orthogonal 
sum of operators defined on finite dimensional (orthogonal) subspaces of 

2.) Moreover, using the matrix of H]^ the matrices of the suboperators can 
be determined explicitly. 

If this proposition is verified, we have obtained a block-diagonalization of 
H^. Thus a way is opened, depending on the numbers a and N, to calculate 
the eigenvalues of with the help of numerical methods. Purely analytical 
solutions of the eigenvalue problem of are also possible, if a = 2, 3, 4. 
But, in these cases one cannot expect that H]^ is a good approximation of 
a realistic Hamiltonian H^. 

6.1.2: In order to verify Proposition 6.1 some further notation is used, which 
is provided by 

Definition 6.2: 1.) Let be given an n G BZF^ and a natural number a, 
which for the moment is completely arbitrary. Then 

(6.1) (n,a) := (ni, • • • ,na) and (a, n) := (n^+i, nQ,+2, • • • )• 

For the infinite second part of h also the abbreviation [a,n) =: f is used. 
2.) Let f be given. Then Nj3{r) denotes the set of all h G BZFn, for which 
(a, n) = f, X^^Q,_|_i rig = N — 13 and < /3 < min{a, N}. Hence the finite 
sequence (n, a) contains exactly f3 numbers 1 and a — /3 numbers 0. Because 
/3 is determined by f, the notation is a bit redundant, but it turns out to be 
useful. 

Now it is supposed that a and N are fixed numbers. Then, Definition 6.2 
yields the following 

Consequence 6.3: 1.) The set Nii{f)^ < /3 < mm{a, A^} is finite, more 
precisely, card ^^{f) = (p. Therefore the subspace T-L^{f) of spanned 
by the ^~^{n) for n G ■N'jiif') has dimension (^). 

2. ) The sets Mp{r) and Npi{r') are disjoint if r 7^ f'. Moreover (3 = /?', if 
and only if the sequences f, r' contain the same number of elements 1. Then 
for each h G J^pir) there is an h' G M/3{'r') such that {h',a) = (n, a). 

3. ) Let B]^ be the ONB of defined by (|A.1.13p . and let ^-^(71) G B^. 
Then there is exactly one (3 so that h G Np{a^ h). Hence the sets Np{f) with 
f = (n^+i, ?T-a+2) ■ ■ ■ ) containing N — (3 numbers 1 and < f3 < min{a, N} 



24 A NEW APPROACH TO THE N-PARTICLE PROBLEM IN QM 



form a complete disjoint dissection of the set of all h £ BZF for which 
i:,n, = N. 

4.) From the above parts 2 and 3 one concludes that the spaces ^^{f) are 
orthogonal for different f, and that they span l-L^, i.e. 

(6.2) n^ = ^n^{f). 

f 

Later on a restriction of the set D of all difference sequences (cf. Notation 
3.8) is needed. 

Definition 6.4: Pq, is the set of all d G D, for which = g > a. In 
addition let Vja := Vj Ci Pa, j = 0, 1, 2. 

Therefore the sets T>ja are again disjoint, and Poa = = {o}- 
6.1.3: Finally a lemma is proved which is basic for the further considera- 
tions. 

Proposition 6.5: If fi G A/'/3(r),0 < /? < min{a,N} and if d G T)a,d / 6, 
then either n + d E J^pir) oi n + d ^ BZF. 

Proof: The proof is complete if one can show that n + d S A/'/3(r) is equiv- 
alent to n d G BZF. First, if n d E A/'^(r), then h + d e BZF holds. 
Second, it follows from h + d ^ BZF, that the (two or one) numbers 1 in 
d must be at positions where there are in h. Likewise, the (two or one) 
numbers -1 in d must be at positions, where numbers 1 are in n. Because 
of d S Va, the sequence (n -|- d, a) has the same quantity /3 of numbers 1 as 
(n, a) has, and f := (a, h) = {a,h + d) because d does not affect f. Thus 
n + d G Mp{f), so that the proof is complete. 

6.2 General properties of the matrix of HJ^ 

6.2.1: In what follows the definitions and results of Section 6.1 are applied 
for the special choice a = a with a properly chosen a. Moreover, the ONBs 
B2 and Bjj are those which are defined via the Hartree-Fock procedure in 
the Sections 4.2 and 4.3. 

6.2.2: In this subsection the first part of Proposition 6.1 is proved. In order 
to do so, the matrix representation of Hjj with respect to is used. The 
proof is complete, if one shows that for h G J^jiir) and h' G A/'/3'(r') with 
f 7^ r': 

(6.3) {h'\H]:j\h) = 0. 

There are three possibilities for the pair n', n. 

If n' — n ^ D, it follows from Formula (14. 3p that (16. 3p holds. 

If n' - n G P\Pq, by Formula ([Oil it is seen that (fOj) holds, too. 

If n' — n G Va, then n' = n -|- d G BZF , d G and d 0. Thus it 
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follows from Proposition 6.5 that h' £ Np{f). This result contradicts the 
supposition h' G Mj3i{'r') and r' 7^ f. Therefore h' — n ^ Da is true. 
Hence ()6.3p holds if {a,h') 7^ {a,h). 

If one denotes the restriction of to the space 'H^{f) by H^{r) one 
obtains 

(6.4) i7^ = 0i?^(f). 

f 

Thus, part one of Proposition 6.1 has been proved. The proof of the second 
part is postponed to Section 6.3. 

6.2.3: In this subsection therefor some preparatory work will be done. Since 
the dimension of 'H^{f) is (^) with < /3 < min{a,N^, it can become gi- 
gantic depending on a and j3. Therefore, aiming at the diagonalization of 
H^{f) it is of vital interest to know how many matrix elements {h'\Hjj{f)\n) 
vanish on principle grounds. This means, how many matrix elements of HJ^ 
are zero for arbitrary dummy Hamiltonians -ff^ respectively their matrices 
E. Then in addition further matrix elements of HJj can be zero for special 
H2Q- But the last aspect will not be considered in this paper. 

Now, from Consequence 4.3 one immediately draws 
Conclusion 6.6: Let n',n S M/3{r) and let n' — h ^ T>a, then 

(6.5) {h'\H-\n) = 0. 
Since by Proposition 3.9 the relation 

h' = h + di + ■ ■ ■ + di, dj eVi 

holds, one has a simple criterion to decide whether n' — h ^ T>a or not. 
Especially, if L > 2 Formula ()6.5p is true. 

6.2.4: Finally the nondiagonal matrix elements with h' — h = d £ Da, d ^ 
are considered. For this purpose let us introduce the following 
Notation 6.7: If n',h G BZE^. Then n',n are called Pjq— concatenated, 
J = 0, 1, 2, if there is a d G T^ja so that h' = h + d. The hzf n', h are simply 
called Da— concatenated if they are Dja— concatenated for j = or 1 or 2. 
With the help of this notation we arrive at the 
Result 6.8: 1.) For each n £ N(j{f), there are exactly 



Ti(a,/3) :=/3(a-/3) 

Piq— concatenated h' G M/s^r). This is because each number 1 out of the /3 
numbers 1 in (n, a) can be transposed at each position of the a — 13 numbers 
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by a d G Pia- 

2.) For each h G ■N'p{f) there are exactly 



2 J \2 

'D2Q-concatcnated h' € Mpif), where (2) = if = 0, 1. This holds because 
each pair of numbers 1 out of the /3 numbers 1 in (n, a) can be brought at 
the position of each pair out of the a — (5 numbers by a d G T>2a- 

3. ) n',n G Mjsif) arc "Doa" concatenated exactly if n' = h. 

4. ) For each h G ■M^{f) there are exactly 

(6.6) r(a, /3) := r2(a, /3) + ri(a, /3) = f ^1 « - f ^1 + - /3) 



Pa-concatenated n' ^ h. 

5.) Finally, let us consider the matrix of H^{r) with the elements {n'\H^\h), 
h',n G Mffif). Then, in the n'-row there are exactly r(a,/3) nondiagonal 
elements which can be unequal zero, and similarly for n-columns. Conse- 
quently, the number Z{oi,(3) of zero nondiagonal elements in each n'-row or 
n-column is 

(6.7) Z{a,^)=l^^ -r{a,^)-l. 

According to Conclusion 6.6 Z{a, /3) is the number of hzf h for which 
h' — h ^ Va in any n'-row, and likewise for the n-columns. 



6.3 The matrices of the operators H^{f) 

In this section the second part of Proposition 6.1 will be proved. This runs 
as follows. 

6.3.1: (3 = a < N. The number z of elements n G Naif) is (^) = 1, and 
the element n has the form 

(6.8) n = (1, • • • , l,n<5+i, • • • )• 
Consequently the matrix of Hjj{f) is of order one and its element is 

(6.9) {h\Hj}\h) = S{n, o) = S{h, o). 

6.3.2: 13 = a — 1 < N. 1.) The number z of elements h G Ma-i{f) is 
(-"^) = a, and the (n, a) for n G A/'a-i(r) contain only one and a — 1 
numbers 1. 

2. ) The n G Na-i{r) are numbered by n =: n^ if is at position a — k in 
(n, a), and k = 0, • • • , a — 1. 

3. ) Any two elements n',n G A/'a-i(^) with n' 7^ n are PiQ-concatenated. 
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This is because the number of n G Ma~i{r), which are 2?a-concatenated 
with h' 7^ h, according to (j6.6|) is 

(6.10) T{a, a — 1) = Ti{a, a — 1) = a — 1 = z — 1. 

4.) The matrix of H^{f) in the present case is a z x z = a x a matrix, 
which has the elements 

(6.11) {k'\H^\k) := {n^i\H^\hf,) = £{h^,hi^' - h^). 

6.3.3: P = a — 2 < N . 1.) The number z of elements n G Na-2{'f') is 
(-"2) = \6L{a — 1), and the (n, a) for n G J\fa-2{'r) contain two numbers 
and a — 2 numbers 1. 

2. ) The n G A/'a_2(^) are numbered by n =: n^x, k < A if the two are at 
the positions a — A and a — n, 0<n<X<a — 1. 

3. ) Any two elements h',h G A/'a-2(^) with h' ^ h are Po-concatenated. 
This is because the number of n G Ma-2{'f'), which are Pa-concatenated 
with h' / n, according to (j6.6|) is 

(6.12) r(a,a - 2) = - 1 = z - 1. 

4. ) The matrix of H^{f ) then is a z x z matrix with z = ^a{a — 1), which 
has the elements 

(6.13) {k',X'\H-\k,X) := {n^'x'\H-\n,^x) = -^("-ka^k'A' -^^a)- 

6.3.4: (3 = 2 < N . 1.) The number z of elements n G M2{r) is (2) = 
^a{a — 1), and the (n, a) for n G A/'2(r) contain two numbers 1 and a — 2 
numbers 0. 

2. ) The h G M2{r) are numbered by h =: fii^x, k < A, if the two 1 are at 
positions k and A, 1<k<A<q; 

3. ) Any two n',n G A/'2(f) are Pa-concatenated. This follow via the same 
argument as in 6.3.3. 

4. ) Like in 6.3.3 one obtains the z x z matrix of H^{'r). It has the elements 

(6.14) {k'\'\H-\k\) := {n^'x'\H~\h^,x) = -^(^^A, "k'A' -n^x)- 

6.3.5: (3 = 1, N > 2. 1.) The number z of elements n G A/i(r) is (") = a, 
and the (n, a) for n G A/i(r) contain one number 1 and a — 1 numbers 0. 

2. ) The h G A/i(f) are numbered by fi =: h^, if 1 is at position k, 1 < k < a. 

3. ) Any two h', n G A/i(f) are PiQ-concatenated. The argument is the same 
as in 6.3.2. 

4. ) Also as in 6.3.2 the matrix H^{f) is obtained. ItisazX2; = axa 
matrix having the elements 

(6.15) {k'\H^\k) := {n^>\H]y\hi^) = i{hi^,h'^ - n^). 
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6.3.6: /3 = 0, > 2. The number z of elements n G 7Vo(f) is = 1, and 
the element h has the form 

(6.16) n = (0, • • • ,0,na+i, • • • ). 
Consequently the matrix of H^{f) is of order 1 and its element is 

(6.17) {h\Hjj\h) = £{n^d) = £{n^o). 

6.3.7: 2 < /3 < a - 2, /3 < iV. 1.) The number z of elements h ^ (r) 
is (^) . It is larger than the numbers z in the previous cases. Each element 
n G ■N'p{r) contains in (n, a) at least three numbers 1 and three numbers 0. 

2. ) For each j3 not all pairs n',n' G -^(^{r) are Po-concatenated. To prove 
this proposition it suffices to give an example. Thus, let 

(6.18) 

h' = (1, 1, 1, • • • , 0, 0, 0, 77.(5+1, • • • )' n = (0, 0, 0, • • • ,1,1,1, Ua+i ■■■), 
and let dj G P, j = 1, 2, 3 be defined by djg = 6jg — 6a+j~3,g, ^> G N. Then 

(6.19) n' = h + di + d2 + d3 

so that h' — h ^ Da- The factual number of non concatinated elements can 
be calculated from Z{a,/3) as defined by Formula (16. 7j) . 

3. ) The elements h G Mp{r) are numbered by n = n^, k = 1, • • • ,z arbitrar- 
ily. Then the matrix elements of H'j^{f) in the present case are 

{k'\H]^\k) : = {h^i\H^\h^) 

(6.20) = i{h^,h^' - hf.) ^n^'-n^GVa 

= , n^/ - n«; ^ Pa . 

6.3.8: Besides the above properties of the matrices of H^{f) the following 
result is of practical relevance. 

Proposition 6.9: Let the sequences r and r' have the same number N —(3 of 
elements 1. If ni,h2 G A/"^ (f) and hi ^ h2, there are h'i,h2 G (f'),n'^ ^ 
h'2 such that 

(6.21) (n;|i?^(f')|n'2) = (ni|i?^(f)|n2), 

and vice versa. Thus, the matrices of H^{f') and H^{f') have the same 
nondiagonal elements. 

Proof: For given rii,n2 the bzf h'l^h^ are chosen according to Consequence 
6.3 as follows: (n^-, a) = {hj, a),j = 1,2. Thus, n'g — h'l = h2 — hi. 
If h2 — hi Pa, it follows from the proof in Subsection 6.2.2. that Formula 
(j6.2ip holds, because both sides are zero. 

Now let us assume that d := h2 — hi G Pq. Then applying Formula (j3.63p 
yields 

(6.22) {hi\H^{f)\h2) = {hi\H^\h2) = C(ra2, m + d, m)E{rh + d, rh), 

rh 
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where the sum runs over ah m with h2 — rh €z BZF and k := m + d ^ 
BZF2. Because E{rh + d, m) = 0, if A:^ = 1, A > a or = !,«:> a, in 
Formula (j6.22p only such components fcp, m^j are relevant, for which p,a < a. 
Therefore, if the condition n2 — rh £ BZF is satisfied, then also — m G 
BZF holds. The other condition is also satisfied, because d = h'2 — h'l. 
Finally, the inversions, which determine C{h,rh + d,rh), only refer to the 
elements of (7^2, a) = (n2,a). Thus one obtains 

(6.23) C{h,m + d,m) = C{h' ,m + d,7h). 
Hence the last term in (j6.22p is equal to 

(6.24) Y,C{n'2,m + d,7h)E{7h + d,m),= {h[\H-\n'2) = (n'i|iJ-(fO|n^). 

m 

The proof of the inverse runs the same way. 

6.4 Conclusion 

The decomposition of the matrix of into orthogonal finite matrices as 
described in the Sections 6.2 and 6.3 now allows, depending on N and a, to 
determine parts of the spectrum of H]^. Thus, it is only a question of the 
capacity of the computers available, which parts of the spectrum one can 
calculate, and it is a question of physical relevance, which parts one wants 
to calculate. Intuitively, the case /3 = a is a Hartree-Fock approximation. 
This suggests that better approximations are achieved if a is greater than 
N, because then the case /3 = a cannot occur. 

7 Final Remarks 

7.1 Summary of the results 

7.1.1: The way of approaching the eigenvalue problem of the Hamiltonian 
presented in the Chapters 3 to 6 may at the first glance seem compli- 
cated. Therefore it is useful to realize the simple kernel of the procedure. I 
will present it in the form of a work program which comprises six steps. 
1** step: As a starting point one formulates the Hamiltonian Hi\f to be con- 
sidered. This is usually done making use of the position-spin representation, 
i.e. Hn is an operator in the Hilbert space = l^^{L^{R^) ®5^). The 
most general form of Hn for charged particles is due to Breit. It can be 
found in the literature, e.g. in [14] p. 247. 

2"^^ step: According to Formula ()2.5p one shapes the dummy Hamiltonian 
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H2Q := -^2^(70), 7o^^ = — 1, belonging to H^. Then one determines the 
orhonormal system Oi via the Hartree-Fock procedure for -ff^, an equiva- 
lent method or an approximation. Eventually, the operators H]^ and H2Q 
are restricted to the spaces H^jH^ with = span 0\ and are denoted 
H~ and H^q. 

S^'^ step: In order to determine the matrix E of one has to choose an 
ONB Bi of TH}. As explained in Section 4.3, the best choise is Bi = Oi. 
Then, using the ONB B2 C 71^ the matrix elements E{k^ m) are calculated 
(cf. e.g. (|A.l.in[) . (|A.1.13|l l. 

4*^ step: From the matrix E one obtains the truncated matrix E by re- 
placing the "small" elements E(k,'m),k ^ rh of E hy zeros as described 
in Section 4.4. The matrix E depends on a number a and determines the 
operator i?^. 

5*^ step: The matrix elements £{n,d) of H^ = r2^(i?2o) calculated 
from the matrix elements E(k,'m) the same way as the £{h,d) are calcu- 
lated from E(k, m) in Section 3.5. In this connection the general results of 
the Sections 6.1 and 6.2 are useful. 

6*^ step: One determines the orthogonal submatrices of the matrix of Hjj 
according to Section 6.3 and diagonalizes them as many as possible numer- 
ically or analytically. 

Then one can try to obtain error estimates applying suitable results of the 
theory of spectral approximation. 

7.1.2: As mentioned at the end of Section 6.4, the lowest energy levels of 
Hjj are not simply Hartree-Fock-like approximations of the true values, if 
is smaller than a. Thus in this case the method presented here reproduces 
also those results, which are obtained by other methods like density func- 
tional theory (DFT) or configuration interaction method (CI). Summing up, 
it is intuitively clear that the approximation of Hjj by H^ is the better the 
smaller the number A^ and the larger the parameter a, which in turn is 
limited by the capacity of computers. My colleague Arno Schindlmayr is 
preparing an application of the proposed method. 



7.2 Finite procedures 

7.2.1: The program described in Section 7.1.1 is a work program in a strict 
sense only if all its steps could be carried through in finite time. Thus, the 
critical points are found in those steps which contain infinite tasks. The first 
and decisive one is the determination of the infinite ONB Oi in the second 
step. For one has to expect that in most cases Oi can neither analytically 
nor numerically be calculated completely. Hence, what can be performed 
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is the determination of a finite part Ofi of Oi, i.e. its elements up to a 
number R. 

However, if only Oji is available, the method described here does not 
brake down. Rather the work program formulated in Section 7.1.1 can also 
be carried through for Ofi instead of Oi. The only question is, which of the 
obtained results are of physical interest. 

In order to get an answer let us use the following obvious notation: 

Moreover, as already introduced above, R is the number of elements of Bfi. 

Now, because the reduced work program depends on the three parameters 
N, R and a, the above question can be answered as follows. 

1. ) If > i?, there are no vectors unequal zero in 'H^_. Therefore this case 
has to be excluded. 

2. ) If N = R, the Hilbert space 'HJ_ is 1-dimensional, so that the A'^- 
particle Hamiltonian HJj^ = HJj^ has one eigenvalue of Hartree-Fock type. 
This result is of minor interest. 

Thus, in order to get better results one needs an R which is "sufficiently" 
larger than N. 

3. ) Let a < N < R. Then, according to the classification described in 
Section 6.3, for each (3 satisfying max{Q, a + N — R}</3<a at least one 
submatrix of the matrix of Hjj^^ exists, which in principle can be diagonalized 
numerically. 

4. ) Let N < a < R. Then, as in point 3, for each /? with max{0, a+N—R} < 
(5 < N there is again at least one submatrix of the matrix of Hjj^, which in 
principle can be diagonalized numerically. 

5. ) If < a = R, then /3 = N. Therefore, there is only one submatrix of 
the matrix of Hjj^, which is identical with the matrix of Hjpj- Moreover 
Hj]\f = Hjj^. Hence, this case is optimal, but it possibly can not be treated 
numerically because R is too large. 

The result of the above considerations now reads: only the cases 3., 4- 
and 5. can be of physical interest. 

7.2.2: Thus, the question arises, are they. In other words, what can be said 
about the spectrum of by studying Hjj^. The answer is given by 

Proposition 7.1: The spectrum of HJ^^ is contained in the spectrum of 

H^. Thus, the finite work program does not change the eigenvalues of H^, 
rather it delivers only a subset of them. 
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The proof runs as follows. It suffices to show that the matrices of Hjj^ and 
H]^ with respect to the ONB Bjj^ are identical. 

Let AifL he the set of all / € BZFl such that = 0, if > R. Then 
fi £ M.fN exactly if ^~^{n) G ^^n- Now let k,7h £ BZF2 and h',h £ M. jn- 
If n — m G BZF and h' — k G BZF, then k,m G Aif2- This implies 

(7.1) Ef{k,m) = E{k,m). 

In addition let T>f be the set of all d := k — ifi, k,ifi G Aif2- Thus dg = 0, 
if ^ > -R, for each d G Vj. 

By definition, £f{h,d), where n G AijN and d G Vj, is constructed 
from Ef{k, rh) via Formula (j3.63p like f (n, d) from -^(A;, rh) (or iS(n, d) from 
E(k,7Ji)). Therefore, by Formula (j7.ip one obtains 

(7.2) £f{n,d) = £{h,d) 
so that by use of (|3.6ip 

(7.3) {n\Hjj^\n) = {n'\H]^\n). 

Finally, if n G AijN and n G J\fp{f), f = {a,h), then Mp{f) C AifN- 
Therefore, corresponding submatrices of Hjjq and Hjj have the same shape. 
Hence, they are identical. 

This result guarantees the practical applicability of the finite work program. 

A Appendix 
A.l Glossary 

A. 1.1 The formalism briefly described in this section was mainly developed 
by Cook [lOj and by Schroeck [llj . The purpose of this appendix is fixing 
notation and formulating some few results, which are used throughout the 
paper. 

The starting point is an axiom of QAI that reads: Let TH} be the Hilbert 
space of a system containing only one particle of a certain kind. Then 
the Hilbert space of a system containing particles of the same kind is 
the symmetric or the antisymmetric subspace of the A^— fold tensor product 
._ (g)^-^i. Likewise the Hilbert spaces of systems composed of different 
kinds of particles are subspaces of appropriate tensor products of one-particle 
Hilbert spaces. 
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In what follows the tensor product (g) of Hilbert spaces is understood to be 
a complete space. But the noncomplete tensor product of linear manifolds 
is a noncomplete linear manifold. This product is denoted ®. 

The inner product in , denoted ( •, • ) or ( •, • )Ar, is defined as usual 
by the inner product ( •, • )i in in the following way. If / = /i (8) • • • (8) /at 
and g = gi<^ ■ ■ ■ <^ gN, then 

(A.1.1) {f,g)N = • • • {fN,gN)i ■ 

By linear and continuous extension {■,-)n is defined onV.^. 

The tensor structure of the iV-particlc Hilbert spaces implies the following 
Proposition A.1.1: Let Bi := {4>x : A G N} be an orthonormal basis 
(ONB) mH^. Then 

(A.1.2) Bn := {(t>M (8 • • • (8 <^A^ : A,- G N,i = 1, . . . , AT} 

is an ONB in . 

Throughout this paper the abbreviation is used: 
(A.1.3) <^Ai-A^ = '^Ai<8)---<8<^A^. 

A.1.2: Let Sn be the symmetric group, and let P G iSjv- Then the operator 
U{P) of the exchange of particles is defined by 

(A.1.4) U (P)<^«i...«^ = <^«^_i(^,...«^_i(^j . 

and by continuous linear extension. 

The operator U (P) has the following properties. 

Proposition A.1.2: 1.) U{P) is invariant under a change of the ONB. 
2.) U{P) is defined on and is unitary, i.e. U{P)U*{P) = 1. Moreover 

(A.1.5) U*{P) = U{P-^), U{PQ) = U{P)U{Q). 

With the help of U{P),P G Sn the symmetrizer and the antisymmetrizer 
are defined by 

(A.1.6) ^^ = W.^ <^^{P)U{P) 

■ P^Sn 

with a~^{P) = 1 and cf~{P) = (— l)"''^'^^ , where J{P) is either the number of 
inversions of P or equivalently the number of transpositions forming P. 

Some useful properties of the operators are summarized in the next 
Proposition A.1.3: are projections defined on V,^ . Moreover 

U{P)S% = S%U{P) = a^{P)S% 



(A.1.7) 
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Then the physicahy relevant subspaces of H are = S^[n^] where + 
stands for bosons and — for fermions. Thus 



(A.1.8) = n1®n^®n^, 

where © is the orthogonal sum as usual. 

A special role in this paper play some orthonormal bases of H^, which 
are defined by 

Proposition A. 1.4: 1.) The set of all vectors 

(A.1.9) := -==5+<A.....«. 

V J '^j ' 

with Ki < . . . < kn and = T1!1=i^k. k.^ is an ONB in H^. Moreover 
2.) The set Bj^ of all vectors 

(A.1.10) := ^/iV!5^,/.,,...,^ 

with Ki < . . . < ktv is an ONB in l-L^. 

A. 1.3: For the problems to be treated in this paper notation ()A.1.9P and 
(jA.l.lOp is not optimal, it can be improved by introducing the sequences of 
occupation numbers by the following 

Definition A. 1.5: 1.) Let be given a sequence of indices , ktv as 

in (lA.l.Qj) or in (lA.l.lOp . where Kj G N,j = I,-- - ,iV. Then define the 
occupation number of k G N by 

N 

(A.1.11) n^ = ^5^^., 

and the sequence of all n^, k £ N, abbreviated 62:/, by 
(A.1.12) (ni,n2,---) =:n. 

2.) Moreover let us denote the set of all hzf ^ for + or for — , by BZF. Then 
the proposition "n is a sequence of occupation numbers" is abbreviated by 
"n is a bzf" or by "n G BZF"^. Sometimes it is useful to write BZFl for 
the set of all bzf I with Ig = L. 

Consequence A. 1.6: 1.) Each sequence of indices ki, • • • , kat determines 
uniquely a hzf , and vice versa. 

2.) The elements of the ONB B% can be written this way: 
(A.1.13) =■■ 
This notation turns out to be very useful. 

A. 1.4: In the next step the question is to be answered which are the physi- 
cally relevant operators, i.e. the relevant observables in 'H^ . Obviously only 
those are relevant which leave the spaces invariant. 



r 1 
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Therefore we define: Let Da C H be the domain of a selfadjoint operator 
A. Then, if for each / S VaI^'H^ the relation Af G holds, the operator 
A is called a physically relevant observable. 
Consequence A. 1.7: 1.) A is physically relevant, exactly if 



(A.1.14) A = S+AS+ + S^AS]^ + Sl^AS'j^ with S*" = 1 - 5+ - 5^ 



holds, i.e. if A is invariant under permutations of particles. 
A. 1.5: Many relevant physical observables are defined using tensor products 
of operators in Hilbert spaces. In order to avoid unnecessary complications 
here only the tensor product of two operators is introduced, because the 
extension to more than two factors is straightforward. 

Thus let be given two Hilbert spaces Tii and and two densely defined 
closed operators Ai and A2. Then A^ and A2 exist having domains P^* 
and Da* ■ 

Now form the (noncomplete) tensor product Vq := Vai'^Pa^ C %i®%2- 
It is dense in 1-Li ® %2 and contains only finite linear combinations of the 
form 



A very useful tool is given by 
Proposition A. 1.8: For bounded operators Ai,A2 the above definition of 
Ai (g) A2 is equivalent to 



where {cpl : X £ N} is an ONE in Tig, g = 1,2 and g = axnfpl 

ni0n2. 

A. 1.6: The results of the last section now are applied to transfer observables 
of M-particle systems into observables of A^-particle systems, M < N. Some 



2.) A is physically relevant, if for each P £ Sn 
(A.1.15) A = U{P)AU*{P) 




Finally the tensor product of Ai and A2 is defined by 
(A.1.18) A^(E)A2=To{At,A*2r. 



Hence, for selfadjoint operators definition (jA.l.lSp reads 

Ai^A2=To{Ai,A2y. 



00 



(A.1.19) 
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basic results in this connection are contained in the fohowing 
Proposition A. 1.9: 1.) Let Am be selfadjoint in H^'^ and let 1 be the 
identity operator in ?^^~*^.Then 

(A.1.20) {Am ® 1)'' = Am®1. 

2. ) If A^j is bounded, then 

(A.1.21) II Am 1 11 = 11 Am \\ ■ 

3. ) If (Am + BmY = All + Bli^ then 

(A.1.22) {Am + Bm) ^ 1 D {Am + {Bm ® 1) • 

The =-sign holds if the domains of both sides are equal, which is the case if 
Am,Bm are bounded and have domain T-L^^ . 

4. ) If Am or Bm is bounded then 

(A.1.23) {A^^1){Bm^I) = {AmBm^I) ■ 

5. ) From ()A.1.7P one concludes that 

(A.1.24) C5 1) = 5± = (5± C5 1)S^ . 

A. 1.7: The operator which defines the physically relevant transfer from T-L^ 
to n'^, M < A is given by 

(A.1.25) nN{AM) := {M\{N - M)\)-^ ^ U{P){Am 01^ ■ ■ .0l)U*{P), 

PeSN 

where Am is a densely defined closed linear operator inT-L^ and where 1 is 
the identity operator on T-L^. It has the following properties. 
Proposition A. 1.10: 1.) The equation 

(A.1.26) U{Q)nN{AM)U*{Q) = J^^(Am) 

holds for each Q £ Sj\f so that Qn{Am) is indeed physically relevant if it is 
selfadjoint. 

2. ) If Am is bounded and selfadjoint with domain the operator ^lj\f{AM) 
is bounded and selfadjoint with domain T-L^ . 

3. ) If Am is unbounded and selfadjoint then Oat (Am) is not necessarily 
selfadjoint. But it is symmetric if it is densely defined. 

A. 1.8: Though the operators ^^^{Am) are physically relevant, they are 
not of interest in a strict sense, if one wants to consider only systems with 
one kind of particles as is the case in this paper. Then only the spaces 
are of interest and the operators defined therein. Thus the 
following definition is natural for the operators of proper physical relevance: 

(A.1.27) I^^(Am) = S^Qn{Am)S^ . 

As in the preceding sections here 0,^{Am) is studied only for selfadjoint 
Am- But the definition itself is much more general. The operators Q^{Am) 
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have some properties, which are of special interest in this paper. 
Proposition A. 1.11: 1.) For each P G Sj\f the following relations hold: 

n%iAM) =Qs^{Amc^i---(e>i)s^ 

(A.1.28) 

2. ) r2^(^jv/) is selfadjoint, because Am is selfadjoint by supposition and 
because are projections, i.e. are bounded. 

3. ) If Am is bounded, then ^^'^{Am) is bounded and 



M 



(A.1.29) ||< (^) \\A 

4. ) If Am + Bm, Am and Bm are selfadjoint. Formula ()A.1.22p implies the 
relation 

(A.1.30) n^A^M + Bm) 3 n%{AM) + n%{BM) ■ 

If Am, Bm are bounded with domain the = sign holds. 

5. ) Let Am+Bm andn%{AM)+^'^{BM)heselfsidiomt. Then (jAXHO]) is an 
equation, because a selfadjoint operator cannot have a selfadjoint extension. 

6. )Finally, from (|A. 1.241) and (|A.1.28p one concludes that 

(A.1.31) n%{n^,{AM)) = n%{AM) . 



A. 2 Dummy Hamiltonians 



A. 2.1: In order to formulate explicitly the dummy Hamiltonians for elec- 
tronic systems it is advisably to work with representations of Hilbert spaces 
instead of the abstract versions used elsewhere in this paper. For our pur- 
poses the position-spin representation is most useful. Therefore the Hilbert 
spaces we are working with are 

N 2 

(A.2.1) cn^ = i^{l\r^)^s^), nlcn^ = i^{l\r^)^s^), 

where is the complex vector space of spin functions u : {1,-1} — )• C 
which is spanned by the ONE (5_i^s}.(Cf. also Subsection 3.1.1.) This 
choice fits into the abstract formulations by the following definition of the 
tensor product. Let x E M, s G {1,-1}, and z := (x,s). Then, if /i,/2 G 
= L2(r3) 51^ one defines /i /a by 

(A.2.2) (/i f2){zi,Z2) = fi{zi)f{z2) . 

For the sake of simplicity let us assume that with the following two examples 
the external fields and the interactions are electrostatic. This means that 
all influences of magnetism and spin are disregarded. 
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A2.2: On the above assumptions the Hamiltonian for the N electrons in an 
atom reads 

N N N 

^ j=l j=l 3 j^k ^'^ 

where Vj = \xj\ and rjk = \xj — x^]- The domain of H^ is as defined 
in Section 3.1.1. Consequently, the atomic dummy Hamiltonian describes 
"dummy helium" and is defined onTi^. It reads explicitly 

(A.2.4) H^, D -^(Pf + P|) - 2joeli^ + ^) + J- 

with 70 = (A^ - l)-i. 

A2.3: Now let us consider N electrons in a finite lattice, the points of 
which are given by y„, a = 1, • • • , AT each carrying the charge cq- Then the 
Hamiltonian is defined by 

N N N N 

j=l j=\ a=l j" j^K 

where r'^^ = \xj — ya\- Thus the dummy solid has a Hamiltonian given by 

and is defined on "H^. 
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